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We present a numerical implementation for virtual corrections to multi-jet production at Next- 
to-Leading order. Using the algorithm of generalised unitarity we compute primitive amplitudes 
from tree-level input. These basic ingredients are then used to compute full colour and helicity 
summed corrections. 
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1. Introduction 

Next-to-leading order (NLO) corrections to observables with multi-particle final states allow 
for precise predictions of the complicated signal and background reactions measured at the Large 
Hadron Collider (LHC) at CERN. Recent years have seen rapid progress in the development of the- 
oretical tools to handle the computation, and automation, of these challenging processes. On-shell 
approaches have been particularly successful at pushing the limit of multiplicity with predictions 
for 2 — > 5 or higher processes being achieved [1-4]. Multi-jet corrections at NLO present an ad- 
ditional level of difficulty in the larger number of parton level processes in both virtual and real 
radiation contributions, nevertheless full NLO predictions have recently been achieved for a four 
jet final state [5]. 

New developments are continuing to improve the algorithms and open up the range and flexi- 
bility of available predictions with a large number under discussion at this conference [6-19] 

In this conference note we present some new developments in the computation of one-loop 
multi-parton processes using the NGluon [20] C++ framework. The numerical library uses on- 
shell methods to compute the basic building blocks of the virtual contributions at NLO which 
we summarise in section ^[ We then look at the construction of full colour and helicity summed 
interference with the born level amplitudes in section || before presenting some results of some 
performance tests in section |J Finally we outline some future directions in the conclusions. 

2. Primitive Amplitudes with NGluon 

The basic building blocks, the primitive amplitudes, of our one-loop amplitudes are computed 
using the generalised unitarity algorithm [21-29] (for recent reviews on the subject see e.g. [30- 
32]) implemented into NGluon . The procedure extracts the coefficients of the scalar integral 
basis from products of one-shell tree-level amplitudes which are computed using Berends-Giele 
recursion relations. Some new features including multiple fermion primitive amplitudes and speed 
improvements through caching and re-using common objects in helicity and permutation sums have 
been described in a recent conference note [33]. The scalar loop integrals functions have been im- 
plemented in a number of public codes [34-36], the default choice in NGluon is FF/QCDLoop. 
The procedure follows in a similar vein to other successful approaches treating multi-gluon ampli- 
tudes [37-41]. 

An «-particle primitive amplitude is defined by an ordered set of n external particles and n 
internal propagators which can be represented by a parent diagram. The full set of independent 
primitives fall into two categories: 

• Amplitudes with a mixture of quark and gluon propagators in the loop, A^f 1 . 

• Amplitudes with an internal fermion loop, A^. 

Each parent diagram can be uniquely specified by a particle in the first position and the propagator 
immediately before hand, a gluon for the mixed case and a fermion for the fermion loop. Each sub- 
sequent propagator will be determined by the external particle and propagator appearing before it, 
either a gluon, quark or blank/dummy propagator in the case that no vertex exists. Some examples 
are shown in figure [[]. 
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The on-shell unitarity method employed in NGluon uses the Four-Dimensional Helicity (FDH) 
scheme. All tree level amplitudes can be computed in four dimensions using the spinor helicity for- 
malism with rational terms extracted from mass shifted fermion and scalar amplitudes. 

A strong cross check on the validity of the rational terms comes from non-trivial relations 
which lead to cancellations in super-symmetric Yang-Mills theories [21, 22, 42]. For the pure 
gluonic amplitudes, or those with a single fermion pair, the rational parts, 2%, obey: 



^(Al ml (l ? ,2,...,(£-1),M£+1)...^ 

-A[ /] (l ? ,2,...,(A:-l),^,(A:+l)...,n)-(-l)"A„ W (l ? ,n,...,(A:+l),^,(A:-l),...,2)) = 0. 



Since the rational terms are the most expensive parts of the computation of the primitive ampli- 
tudes using these relations to reduce the total number of rational terms leads to a considerable 
improvement in evaluation time. 



3. Colour Summed Amplitudes 

The full colour amplitudes can be written as a decomposition of colour factors and partial 
amplitudes, 



where / is the loop order. 

In the above {/} and {a} represent fundamental and adjoint SU (N c ) indices respectively while 
{/?} and {h} represent the momenta and helicity. The one-loop partial amplitudes must be further 
decomposed into the primitive objects described in the previous section. For the pure gluonic case 
and the case of a single fermion pair such a decomposition is known to all-orders. As an explicit 




(2.1) 



(2.2) 



^ W (M, w) = Iq(w, {i})4 } m,{h}), 
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Figure 1: Examples of primitive amplitudes 
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example, the gluonic amplitudes are written as, 



f ] = £ tr(l,a 2 ,...,a n )A^(l : a 2 ,...,a n ) 

oeS„-i/Z 

'P = £ N c tr(l,°2,...,On)A^(l,0 2 ,---,Gn) 

[n/2]-\ 



+ L E tr{\,0 2: ...,C c -i)tr{0 c ,...,0 n )A\r) : {\,02,...,0n), 

c=3 aeS„/S, KC 



(3.2) 



(3.3) 



where 



A^(l,{a},{ j 8}) = (-1)" £ Ag(l,(T 2 ,...,c7 n ) 

creCOPM^} 



(3.4) 
(3.5) 



The sum 'COP' is over all possible mergings of the sets {a} = {2, . . . , c — 1 } and {/3 } = {c, . . . , n) 
while keeping a r fixed and using the cyclically ordered permutations of /3. For the multi-fermion 
amplitudes the decompositions can be constructed by a systematic matching of terms to a Feynman 
diagram representation [31, 43]. 

The NLO corrections come from the interference between tree-level and loop-level amplitudes, 



do 



virtual 



(3-6) 



Since the number of colour structures becomes extremely large with the increasing multiplicity it 
is important to reduce the colour basis as much as possible. We apply photon decoupling identities 
(and their generalisations to the multi-fermion case) in order to reduce the amplitudes to an (n — 2) ! 
of tree amplitudes. As noted by Dixon, Del Duca and Maltoni [44], this has some rather striking 



simplification in the final representation of eq. (3.6). 
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24 
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£/ 6 (d,d,g,g,g,g) 


24 
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59 


£f 6 (d,d,u,u,g,g) 


12 


80 


14 


£/(,(d,d,u,u,s,s) 
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32 
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Table 1: The number of independent primitive amplitudes appearing in the colour sums at tree-level (A„ ) 

[ml Iff 

and at one-loop for the mixed (A„ ) and fermion loop (A„ ) cases. Like-flavour amplitudes for multiple 
fermions are obtained by (anti-)symmetrization and therefore contain larger bases of primitives. 
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4. Performance Tests 



4.1 Infra-Red Poles 



The first important check on the implementation is the verification of the universal Infra-Red 
and Ultra- Violet poles in the dimensional regularisation parameter e [45, 46]. For the interfered 
amplitude in massless QCD these take a rather simple form (in this case unrenormalised in the 
FDH scheme): 

Re^^-FDH^ = _ j_ (N c n g + ^n q ^ K (0) | 2 

+ finite terms, (4.1) 

where n q is the number of external quarks and n g is the number of external gluons. N c is the number 
of colours and «/ is the number of light flavours and j3o = 1 1Nc ^ 2N f _ The quantity |^°^| 2 is the 
colour correlated Born amplitude defined by 

1^(0) f = ^(0),t. ^..^(O). (42) 

where T g = if" bc and T q = T?.,{J q = —TV). Scheme conversion to Conventional Dimensional 
Regularisation (CDR) and renormalisation are all corrections proportional to the Born amplitude 
and are defined as follows: 

Re(^ ^.(^' FDH > R -^ FDH ' U )) _ (n g + n q -2)fp N c 



|^(°)| 2 2 \e 3 

Re(^)* ^)> CDR ' U - ^ FDH > U )) N c n q ( N 2 C 



|^/(°)| 2 3 AN C \ 3 

Re(^°^.(^ cm ' R -^ FDH ' u )) {n g + n q -2)Po N c n q ( N 2 



(4.3) 

1 J (4.4) 



|^(0) |2 2s 3 4N C V 3 7 (4 ' 5) 

4.2 Accuracy 

We examine the accuracy of the amplitudes against the known pole structure and by using a 
re-scaling of the external momenta to test the finite terms. 10 5 phase-space points were generated 
using the RAMBO algorithm using rather weak kinematics cuts requiring only that p t ■ pj > I0~ 4 s 
for a centre of mass energy of s = 7 TeV. 

The thinner lined histograms show the number of points re-evaluated in higher precision 
(quadruple precision using the qd package [47]). In every case, except for the six gluon ampli- 
tude, all re-evaluated points achieved the required relative accuracy of 10~ 4 . The number of points 
for re-evaluation reaches a maximum of 2.2% for six gluons and decreases rapidly with the number 
of fermion pairs. In this most complicated channel 4 from 100,000 events required octuple preci- 
sion to reach the desired threshold, nevertheless they would likely be excluded by common LHC 
cuts. Rough evaluation times for the most complicated 2 — > 4 processes are of order 15 — 20s for 
full helicity and colour and including re-scaling to compute an accuracy estimate (see table 
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Accuracy Accuracy 

Figure 2: Accuracy plots for the 2 — > 4 processes contributing to the 4-jet rate. Relative accuracy of 10 5 
points over a flat phase-space. Points appearing in the shaded region are re-evaluated in quadruple precision. 



n 


4 


5 


6 


w-gluons 


0.03s 


0.5s 


15.0s 


qq+(n — 2)-gluons 


0.03s 


0.75s 


17.0s 



Table 2: Estimated evaluation times (s) for n-gluon and qq + n-gluon amplitudes using Intel Core i7 2.7GHz. 



5. Conclusions 

We have described the implementation of full colour QCD virtual contributions to NLO multi- 
jet production at hadron colliders. The primitive amplitude decomposition and use of on-shell gen- 
eralised unitarity shows to be a fast and accurate method for the evaluation which scales well with 
the number of final state particles. We are looking forward to some phenomenological applications 
in the near future. 
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